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ABSTRACT 

Many clusters of galaxies contain an appreciable amount of hot gas, the intracluster medium. As 
a consequence, gas will be stripped from galaxies that move through the inter cluster medium, if the 
ram pressure exceeds the internal gravitational force. We study the interaction between the intracluster 
medium and an extended gas component of dwarf galaxies confined by a surrounding cold dark matter 
halo analytically and numerically, using axisymmetric two-dimensional hydrodynamical simulations at 
high resolution. The results show that the gas within the dark matter halo is totally stripped in a typical 
galactic cluster. The process of ram pressure stripping therefore must have played an important role 
during the chemo-dynamical evolution of dwarf galaxies in galactic clusters. Our results predict that 
most of the chemical evolution of dwarf galaxies must have occurred at high redshift, before the virialized 
cluster had formed. 

Subject headings: galaxies: dwarf — formation — evolution — clusters — intergalactic medium 
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1. INTRODUCTION 

Explaining the morphology-density relation (Dressier 
1980), that is the higher early-type galaxy fraction in clus- 
ters of galaxies in contrast to the higher late-type galaxy 
fraction in the fields, remains one of the most impor- 
tant problems in cosmology. This relation holds also for 
dwarf galaxies. Gas-poor dwarf galaxies are strongly clus- 
tered, while gas-rich dwarf galaxies appear to be the most 
weakly clustered objects (Binggeli et al. 1987; Ferguson & 
Sandage 1988; Binggeh et al. 1990, Thuan et al. 1991). 

Since many clusters of galaxies contain an apprecia- 
ble amount of hot gas, the intracluster medium (ICM), 
gas will be stripped from galaxies that move through the 
ICM, if the ram pressure exceeds the internal gravitational 
force. There are actually many observational evidences 
for the ram pressure stripping of giant galaxies in clus- 
ters of galaxies (Irwin et al. 1987; White et al. 1991; 
Bohringer et al. 1995). In addition, the physics of ram 
pressure stripping of giant galaxies has been studied in de- 
tail using hydrodynamic simulations (e.g. Balsara, Livio, 
& O'Dea 1994 and references therein). Since dwarf galax- 
ies have smaller escape velocities, ram pressure stripping 
is expected to be more efficient than for giant galaxies. 
However, no attempts have ever been made to examine 
the amount of gas that could be stripped from low-mass 
galaxies with shallow gravitational potential wells during 
their passage through an ICM. 

This situation motivated us to consider the ram pressure 
stripping of the diffuse gas phase of dwarf galaxies by the 
ICM which may be of prime importance to understand the 
morphology-density relation for dwarf galaxies. We esti- 
mate the critical condition for gas removal by the interac- 
tion between an ICM and gaseous dwarf galaxies confined 
by a surrounding cold dark matter (CDM) halo. More- 
over, we verified that condition using an axial symmetric 
two-dimensional hydrodynamic code which is based on the 
piecewise parabolic method (PPM) described by Colella & 
Woodward (1984). 
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In §2, prior to the main subject, we describe the model 
of dwarf galaxies which is treated as a one-parameter fam- 
ily defined by the core mass of the CDM halos, using the 
observed scaling relation (Burkert 1995). In §3, we give the 
analytical estimates of the critical conditions for the gas 
ablation from dwarf galaxies due to ram pressure stripping 
by the ICM. In §4, we demonstrate the complex interaction 
between dwarf galaxies and the ICM using hydrodynam- 
ical simulations and compare the simulation results with 
the analytical predictions. In §5, we summarise the results 
of this paper, and discuss observational and theoretical im- 
plications. 

2. MODEL OF DWARF GALAXIES 

The hierarchical clustering model of galaxies has been 
successful in explaining the clustering pattern of galax- 
ies revealed by redshift surveys. On the other hand, the 
hierarchical model predicts a large number of low-mass 
galaxies formed at high redshifts beyond that estimated 
from the observed luminosity function of galaxies (White 
& Frenk 1991; Cole et al. 1994). Therefore, several mech- 
anisms for suppressing and/or delaying the formation of 
dwarf galaxies have been proposed in order to remove this 
serious discrepancy. 

Using a three-dimensional iV-body/SPH simulation 
code combined with stellar population synthesis, Mori et 
al. (1997) showed that an energy feedback via supernovae 
and stellar winds from massive stars keeps a large fraction 
of the gas hot and suppresses the formation of dwarf galax- 
ies. Moreover, they pointed out that this energy feedback 
is a necessary mechanism to reproduce the internal struc- 
ture and the photometric quantities of the nearby dwarf 
galaxies (see also Dekel & Silk 1986; Yoshii & Arimoto 
1987; Burkert & Ruiz-Lapuernnte 1997; Mori, Yoshu & 
Nomoto 1999). 

Babul & Rees (1992) and Efstathiou (1992) argued that 
the formation epoch of dwarf galaxies is delayed until z < 1 
due to the photoionization of the gas by the ultraviolet 

8588, Japan 
1 



2 



Gas Stripping of Dwarf Galaxies in Clusters of Galaxies 



background radiation at high redshifts. The ionizing back- 
ground at z > 1 is high enough to keep the gas in dwarf 
galaxy halos confined and stable, neither able to escape, 
nor able to collapse (see also Rees 1986; Ikeuchi 1986; 
Katz, Weinberg & Hernquist 1996; Navarro & Steinmetz 
1997). 

Following this scenario, we assume that a (proto-) dwarf 
galaxy has an already virialized CDM halo and a large 
amount of extended hot gas heated either due to a first 
population of supernovae or due to the photoionization by 
ultraviolet background radiation at the formation epoch 
of the galaxy clusters. 

2.1. Structure oj the CDM halos 

The formation of the CDM halos through hierarchical 
clustering predicts that the equilibrium density profile of a 
CDM halo has a central cusp (Dubinski & Carlberg 1991; 
Navarro, Frenk & White 1996; Fukushige & Makino 1997; 
Moore et al. 1998). However, recent observed rotation 
curves of nearby dwarf galaxies rule out the singular den- 
sity profile of dark matter halos (Flores & Primack 1994; 
Moore 1994). Moreover, Burkert (1995) points out that 
the density profile 



[r + ro)(r2 + rg) ' 



(1) 



nicely reproduces the rotation curves of nearby dwarf 
galaxies and the central density pdo is correlated with the 
core radius rg through a simple scaling relation. 

In this paper we assume that the density distribution of 
the CDM halos is represented by Burkert's (1995) profile 
for a wide range of core masses from IQ^Mq to 10^^ Mq. 
For the density distribution of the CDM halo in equation 
(|^),the potential is found by integrating Poisson's equation 
as 
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The central potential is given by $d(0) = —n^GpdorQ, 
where G is the gravitational constant. The mass distribu- 
tion of the CDM halo is given by 
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Using the observed scaling relation derived by Burkert 
(1995), pdo and ro are calculated as 



ro = 3.07 



Mo 



109 M, 



kpc. 



and 
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where the core mass Mq is the total mass of the CDM halo 
inside rg. 

The circular velocity v'^ = GMd{r)/r of the CDM halo 
has a maximum value 



48.7 



Mo 



109 Mp 



km s 



(6) 



at the radius of r^max = 3.24ro (see Fig. 1). This velocity 



and the radius are related by 



10.5 



kpc 



km s 
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2.2. Structure of the gas 



We focus on the general question of the interaction be- 
tween the intracluster medium and a hot gaseous com- 
ponent in dwarf galaxies that is generated e.g. by the 
supernova heating (Dekel & Silk 1986; Burkert & Ruiz- 
Lapuente 1997; Mori et al. 1997; Mori et al. 1999), or 
by the photoionization due to the ultraviolet background 
radiation (Efstathiou 1992; Babul & Rees 1992). In this 
case the gas is pressure supported justifying the assump- 
tion of an initially spherically symmetric distribution and 
of axi-symmetry. 

The self-gravity of the gas is neglected for simplicity and 
the gas is assumed to be in hydrostatic equilibrium with 
the constant temperature 
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where p is the mean molecular weight, rup is the proton 
mass, and fee is the Boltzmann's constant. The density 
distribution of the gas is given by 



{$d(r)-<i>d(0)} 



(10) 



where pgo is the central density of the gas. The mass ratio 
between the gas and the CDM halo within a core radius 
is defined as 



F ■ 



Mo 



(11) 



where Mgo is the total gas mass inside a core radius. The 
relation between the central density of the gas and of the 
dark matter is given by the numerical solution of equation 
(0): Pgo = 1.40i^pdo- We assumed initially F = 0.1 in 
this paper. Consequently each dwarf galaxy is described as 
a one-parameter family by the core mass Mo of the CDM 
halos. 

3. ANALYTICAL MODEL 

The stripping process is generally classified into two dis- 
tinct types such as the instantaneous stripping phase and 
the continuous stripping phase due to Kelvin-Helmholtz 
instability. In this section, we estimate analytically the 
gas ablation affects on the evolution of a dwarf galaxy. 
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3.1. Instantaneous ram pressure stripping 

We assume that the CDM halo associated with the dwarf 
galaxies moves in a homogeneous ICM with number den- 
sity of ncG = Pcg/(m'71p) ~ 10~^ cm~'^ and relative 
velocity Ugai ~ 10'^ km s~^ corresponding to the three- 
dimensional velocity dispersion of the galaxy cluster. 

The complete condition for the ram pressure stripping 
from a dwarf galaxy requires that the ram pressure of the 
ICM exceeds the thermal pressure Pgok-QT / {jjLm-p) at the 
center of the gravitational potential well of the CDM halo. 
Using equation (0), such a condition is described by 



2 ^ GA/oPgo 



(12) 



Thus, if the core mass Mq is smaller than the critical core 
mass 



Mis = 1-27 x 10^ 
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the gas in a dwarf galaxy is totally stripped by the ram 
pressure of the ICM. 

The collision between the gas of the dwarf galaxy 
and the ICM generates a forward shock that propagates 
through the gas of the dwarf galaxy and a reverse shock 
that propagates through the ICM. The time-scale of the 
gas removal is roughly estimated by dividing the diameter 
of the core by the velocity of the forward shock Vfs- The 
analysis of the one-dimensional shock problem gives 



PCG, 

Pso 



-Wgal, 



(14) 



for the high-speed collision of two homogeneous non- 
gravitating media with large density ratio (pcG ^ Pgo)- 
The time-scale of the mass removal from dwarf galaxies is 
estimated by 
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Consequently, the gas is totally stripped from less- 
massive dwarf galaxies {Mq <^ lO^M©) in the typical envi- 
ronment of the ICM because the ram pressure of the ICM 
exceeds the gravitational force of these galaxies. The time- 
scale of the gas removal is small compared with the char- 
acteristic time-scale (several Gyr) a galaxy spends within 
potentials, therefore, lose their gas easily within this dwarf 
galaxies (Mq ^ IO^Mq) the gas might be significantly re- 
moved except around the central region of the gravita- 
tional potential well. 

3.2. Kelvin- Helmholtz instability 



Gas in the massive dwarf galaxies surviving the instan- 
taneous ram pressure stripping is subsequently removed 
by Kelvin-Helmholtz instability occurring at the interface 
between the gas in the dwarf galaxy and the ICM. Mur- 
ray et al.(1993) estimated the characteristic growth time 
of Kelvin-Helmholtz instability for a dense cloud embed- 
ded in a low-density background. They showed that the 
growth time is comparable to the sound crossing time of 
the gas cloud if gravity is negligible. On the contrary, the 
presence of a gravitational field by the surrounding CDM 
halo tends to suppress the instability and to stabilize the 
gas against removal from the dwarf galaxy. 

We suppose that the process of the instanteneous ram 
pressure stripping removes the gas distributed outside a 
core radius tq. Moreover, a sharp discontinuity of the gas 
density is established between the ICM and the galactic 
gas having nearly constant density after the shock passage. 
The unstable wave number of Kelvin-Helmholtz instability 
at the interface for the incompressible fluid is given by 



k> g 



PCG 



,avg 



PCGPg,avgW„al 



(17) 



where 



;.avg 



is the mean gas density inside and g = 



GM() /rQ is the gravitational acceleration at the fluid inter- 
face (cf. Chandrasekhar 1961). Since the dominant wave- 
length for the gas ablation by Kelvin-Helmholtz instability 
is the order of tq (Murray et al. 1993), the inequahty ( pT| ) 
for PCG < 



is transformed as 

2 ^ GAfoPg,avg 



(18) 



Using this inequality, it is found that the diffuse gas com- 
ponent in these galaxies is removed if their core mass Mq 
does not exceed a critical value given by 



Mkh = 1.60 X 10^' ( — 



"CG 



10 *cm^3 

^^gal 



lO^km s" 



M0(19) 



This value corresponds to masses characteristic for a gi- 
ant galaxy which indicates that dwarf galaxies should in 
general be affected by Kelvin-Helmholtz stripping. The 
gas ablation due to Kelvin-Helmholtz instability is effec- 
tive for massive dwarf galaxies if the mass-loss time-scale is 
smaller than the characteristic time-scale a galaxy spends 
spends within a cluster environment. 

According to Nulsen (1982) the mass-loss rate from the 
galaxy through Kelvin-Helmholtz instability is estimated 



as 



Mkh = TrroPcGWgai, 



(20) 



the gas removal time-scale of the dwarf galaxies is there- 
fore given by 



FMn 



(21) 
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Table 1 
Parameters of the I CM 



Model 




^^gal 










cm~^ 


km s"-^ 


K 




erg cm~^ 


(a) 


1.0 X 10-3 


1000 


1.0 X 10^ 


2.1 


1.0 X 10-11 


ib) 


1.0 X 10--'' 


500 


1.0 X 10^ 


1.1 


2.6 X 10-12 


(c) 


1.0 X 10"'' 


1000 


1.0 X 10^ 


2.1 


1.0 X 10-12 


id) 


1.0 X 10"^ 


500 


1.0 X 10^ 


1.1 


2.6 X 10-13 


(e) 


1.0 X 10-5 


1000 


1.0 X 10^ 


2.1 


1.0 X 10-13 



'^Temperature of the ICM 
I'Mach number of the ICM 



Moreover, it should be noted that the characteristic time- 
scale, normalized by the dynamical time 



'''dyn — 



32(j'/3(j,avg 



:8.93xl0^('^^l yr. 



(23) 



(24) 



where Pd,avg is the mean density of the CDM halo inside 
ro, is given by 

TKH 24 ,^ "'CO 



Tdyn 



^"^ (o.l) (lO-^cm-a) 



X I ) (25) 

103km s^i ' ^ ' 



This equation does not have an explicit dependence on 
the mass of the CDM halo. Massive dwarf galaxies lose 
the extended gas within ~ 25rdyn for ncG = IQ-* cm-3 
and Wgai = 1000 km s-i. Consequently, we expect that 
even for the massive dwarf galaxy the extended gas will 
be stripped through Kelvin- Helmholtz instability within a 
short time-scale. 

On the other hand, Balsara, Livio, & O'Dea (1994) 
demonstrated that the gas accretes from downstream into 
the core in their hydrodynamical simulations of the inter- 
action between the ICM and the giant galaxy. Since the 
dwarf galaxies have shallower potential wells than the gi- 
ant galaxies, we neglect the effect of the gas accretion form 
backward into the galaxy. It may play however a role for 
the massive dwarf galaxies where our estimation of the gas 
ablation might not be suitable. 

4. NUMERICAL MODEL 

In addition to the effect of the accretion inflow into 
the core from downstream, the above analytic arguments 
neglect the geometrical effects between the gas in dwarf 
galaxies and the ICM and the presence of the Rayleigh- 
Taylor instability. Moreover, the compressibility of the gas 
may become important because the each galaxy moves in 
the cluster of galaxies with the transonic velocity that is 



corresponding to the velocity dispersion of the cluster of 
galaxies. Therefore, a realistic analysis that uses a hy- 
drodynamic simulation is necessary to examine the effect 
of the gas removal due to the interaction between a dwarf 
galaxy and an ICM. In this section, we describe the numer- 
ical model of ram pressure stripping from dwarf galaxies. 

4.1. Initial conditions and simulation method 

For modeling the interaction between a galaxy and an 
ICM, Portnoy, Pistinner & Shaviv (1993) showed that a 
treatment of the protons and the electrons as two fluids 
resulted in negligible differences for the mass of the gas 
inside the galaxy. In this paper, therefore, the evolution 
of the gas is described by the hydrodynamic equations for 
the single perfect fluid. The continuity equation, the mo- 
mentum equation, and the thermal energy equation are 
given by 



dp 
dt 



+ V • {pv) = 0, 



dpvv 

dt 



V • (pvv) + yp = -v$d 



and 



dpe 
~dt 



V • {pev) + V{Pv) = -pv ■ V$d, 



(26) 



(27) 



(28) 



where p is the gas density, vis the gas velocity, P is the 
gas pressure, 7(= 5/3) is the adiabatic index, and e is the 
total specific energy 



1 2 



1 P 
7 - 1 p ■ 



(29) 



Initial conditions are generated using the descriptions 
in §2. We neglect the effects of the self-gravity of the gas, 
radiative cooling, the heating by supcrnovae and stellar 
winds from massive stars, and the photoionization by ul- 
traviolet radiation background. We discuss them in §5. 
The equations are solved by a finite difference code VH-l. 
VH-1 is based on the piecewise parabolic method (PPM) 
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described by Colella & Woodward (1984) and was writ- 
ten and tested by the numerical astrophysics group at 
the Virginia Institute for Theoretical Astrophysics. Since 
the PPM scheme has a great advantage due to the reduc- 
tion of numerical viscosity, all fluid interfaces are sharply 
preserved and small-scale features can be resolved. This 
scheme is, therefore, suited for this class of problems. 

The system is assumed to be axial symmetric and de- 
scribed by cylindrical geometry. The center of the fixed 
gravitational potential given by equation is located on 
the symmetric axis and the gas distribution of the dwarf 
galaxy is set up using equation . The reflecting bound- 
ary condition is adopted at the bottom boundary that cor- 
responds to the symmetric axis and the ICM flows con- 
tinuously from the left stream boundary parallel to the 
axis. The outflow boundary conditions, by imposing for 
each variable a zero gradient [d/dr = 0), are adopted at 
the top and the right stream boundaries. Unfortunately, 
this condition does not eliminate reflection waves from the 
boundaries. Thus, the simulation box is separated by two 
parts such as the inner rectangle 8ro long and 3ro wide 
and an outer surrounding part to keep the boundaries far 
away. The grids are equally spaced with 100 zones per 
core radius for the higher resolution runs or 50 zones for 
the lower resolution runs in the inner part, and are ex- 
ponentially spaced in the outer part. A useful additional 
effect of this non-uniform grid is the increased numerical 
dissipation of disturbances that propagate to a large dis- 
tance. This helps solving the problem of residual reflection 
waves from the boundary. 

4.2. Results of simulations 

We have performed a parameter study varying the core 
mass Mo of the CDM halo from IO^Mq to 101°Mq, vary- 
ing the relative velocity of the dwarf galaxy Wgai from 500 
km s^^ to 1000 km s^^, and varying the number density 
ncG of the ICM from 10"^ cm'^ to lO'^ cm'^. Table 1 
displays the model parameters of the ICM in this paper. 

Figure 2 shows snapshots of the run for the model (c) 
of Mo = IO^Mq, ncG = 10"'' cm-3, and Vgai = 1000 
km s^^ as a function of elapsed time at 5.61 x 10^,2.86 x 
10^4.58 X 10^,7.42 x 10^ and 9.63 x 10^ yr. The left 
and the right colour images show the logarithmic density 
distribution and the logarithmic pressure distribution, re- 
spectively. Arrows in the right panels indicate the velocity 
vector of the gas flow at each point and their lengths are 
proportional to the absolute values of their velocity. The 
CDM halo is located at the origin and the ICM flows in 
from the left side. 

The upper panel shows the early phase of the interac- 
tion between the ICM and the gas of the dwarf galaxy. 
It is clearly seen that the forward shock that propagates 
through the gas of dwarf galaxy, and the reverse shock 
that propagates through the ICM are formed in the down- 
stream and the upstream of the contact discontinuity re- 
spectively. The propagation of the forward shock causes 
the strong compression of the gas in the dwarf galaxy as 
seen in the upper and the second panel. Below the second 
panel, Kelvin-Helmholtz instabilities are observed at the 
contact discontinuity. 

Since the instabilities are suppressed by the gravity, the 
growing modes are only the disturbances of the large wave 
numbers that have the shortest growth time. Disturbances 



with small wave numbers grow while moving away from 
the galaxy. The third panel shows that the gas associated 
with the dwarf galaxy is accelerated by the ram pressure 
of the ICM and is pushed out of the potential well. Since 
the front of the accelerated gas is highly Ray leigh- Taylor 
unstable, there are many small irregularities in front of 
the contact discontinuity. Accurately, the feature seen 
on the symmetric axis (i? = 0) at ~ 0.5ro is caused 
by pure Rayleigh- Taylor instability. The combination of 
Rayleigh- Taylor instability and Kelvin-Helmholtz instabil- 
ity forms other irregularities because the gas flows tangen- 
tially to the fluid interface as seen in the third and right 
panel. However, the perturbations do not grow signifl- 
cantly because the gravity of the CDM halo suppresses 
Rayleigh- Taylor instability. The gas removal time-scale 
nicely agrees with the analytical estimate in equation (16) 
which predicts tjs = 5.42 x 10'' yr. 

Figure 3 shows snapshots of the run for the model (c) 
of Mo = lOi"M0, ncG = 10"^ cm^^^ and Wgai = 1000 
km s~^ as a function of elapsed time at 6.84 x 10'', 6.99 x 
10^ 1.15 x 10^, 2.18 X 10^, and 3.26 x 10^ yr. 

The early evolution of the interaction between the gas 
and the ICM is almost the same as in the case of Mo = 
1O^M0. However, since the gravitational potential is deep 
and the central thermal pressure is larger than the ram 
pressure of the ICM as shown in the analytical model, the 
mass-loss process is not instantaneous but mild ablation 
due to Kelvin-Helmholtz instabilities occur. Moreover, the 
gas interface is Rayleigh- Taylor stable because the gas ac- 
celeration is week. The velocity field in the right bottom 
panel reveals that the flow turns back and an accretion in- 
flow into the core develops close to the symmetry axis on 
the downstream side. Accretion occurs quasi-periodically 
as a radial-pumping mode. The condition of the flow is 
summarized as follows. 

1. This flow from the downstream side and the post- 
shock gas from the upstream side cause increased 
compression mainly along the symmetry axis. 

2. The gas in the galaxy, therefore, acts by expanding 
sideways as seen in the third panel. Since this ex- 
pansion increases the cross section of the interaction 
between the galactic gas and the ICM, the mass- loss 
rate increases. 

3. The expansion decreases the pressure gradient 
around the galaxy center and the gas contracts to 
the galaxy center by the gravitational force. Then 
the system recovers the quasi-static states again. 

4. This cycle repeats again. 

This process is effective in the simulations of a massive 
dwarf galaxy or of a small ram pressure of the ICM. For 
this specific process, the analytic condition given by §3.2 
fails to describe the gas ablation from the dwarf galaxy. 

Figure 4 shows the evolution of the gas mass inside the 
core radius around the center of the dark matter halo as 
a function of time for core masses of the CDM halo from 
1O^M0 to IO^^Mq. Each curve is categorized by the model 
parameters (see Table 1). The gas is instantaneously 
stripped in low-mass dwarfs or for large ram pressures. 
Due to the accretion inflow into the galactic core from the 
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down stream, the net mass-loss rate becomes small. More- 
over, we can observe oscillations of the total mass inside 
a core radius due to the mass accretion. These oscilla- 
tion modes have periods of several dynamical timescales. 
Even for the massive dwarf galaxy {Mq — lO^^M©), the 
gas could be instantaneously stripped if the ram pressure 
is large as e.g. in the cores of galactic clusters. 

Figure 5 summarizes the results of our simulations. This 
diagram shows the mass-loss rate due to ram pressure 
stripping as a sequence of the core mass. Filled circles 
indicate the cases where the gas in the core radius is com- 
pletely stripped within 1 Gyr. In these cases the ram pres- 
sure stripping by the ICM is very effective and the whole 
gas in the potential well is rapidly removed. Open circles 
indicate the cases where the gas in the core radius is not 
completely stripped within 2 Gyr. In these cases, the ram 
pressure stripping by the ICM is not so effective. The thin 
line indicates the relation of 



2 GMoPgo 

PCGWgal ~ o ' 



(30) 



which is the instantaneous stripping condition described 
in §3.1. This line divides the parameter space into two 
regions. Stripping is effective for the upper region and 
is ineffective for the lower region. This relation roughly 
agrees with the numerical experiments in the range from 
Mo = lO^M© to Mq = 101°Mq. The dashed hue indicates 
the relation of 



27rrn 



(31) 



which is the Kelvin-Helmholtz stripping condition de- 
scribed in §3.2. Stripping by Kelvin-Helmholtz instability 
is effective for the upper region and is not effective for the 
lower region. Since there is mass accretion from the back 
of the galaxy, this relation is not well reproduced by the 
numerical experiments. 

Using the data of intermediate cases that are shown by 
filed triangles in Figure 5, we find a critical core mass for 
effective ram pressure stripping of 



2.52 X 10- 



"CG \ 5 / Wgal 



/ "CG \ 

U0-4cm-3/ 



Vl03k 



m s 



The thick line shows the loci of this critical core mass. 

In Figure 5, rcc is the distance from the center of galaxy 
cluster assuming the isothermal /3-model 



"-CG = '^CGO 



/ rcG 



V'^CGO 



/3 



(33) 



with f3 = 0.6,ncG0 = 2 x 10"-^ cm-^rcco = 0.25 Mpc, 
and acG = 866 km s~^, where ncGo and rcco is the cen- 
tral number density and the core radius of an ICM and 
(TcG is the line-of-sight velocity dispersion of the galaxy 
cluster which corresponds to the relative velocity of the 
dwarf galaxy. If the galaxy number density follows an ap- 
proximate King model as 



"gal OC 



rcGo) 



-3/2 



(34) 



the median radius of the galaxy distribution for r < 
lOrcGo is about 3.5rcG0- The shaded region corresponds 
to ICM conditions inside this radius. The diagram in- 
dicates that the process of the ram pressure stripping is 
very effective for dwarf galaxies in the clusters of galax- 
ies. There, diffuse gas in dwarf galaxies should rapidly be 
removed. 

5. SUMMARY AND DISCUSSION 

The physics of the interaction between the ICM and 
dwarf galaxies confined by a surrounding CDM halo has 
been studied in detail using analytical estimates and hy- 
drodynamic simulations. We have performed a parameter 
study varying the core mass of the CDM, the relative ve- 
locity of the galaxy and the density of the ICM. We find 
that the gas in dwarf galaxies is rapidly removed in a typ- 
ical cluster environment by ram-pressure stripping. 

Our results can be applied to clusters of galaxies, where 
X-ray emission has been observed. Their gas number 
density ncG at the median radius of the galaxy distri- 
bution and the line-of-sight velocity dispersion acG are 
plotted in Figure 6. Data for the central gas density and 
/3 are from Jones & Forman (1999) and Briel, Henry & 
Bohringer (1992), and the velocity dispersions are from 
the same sources including also Hughes (1989). One par- 
ticular example, the Coma cluster with no = 2.89 x 10^"^ 
ciR~^,P = 0.75 and ctcg = 1010 km s~^, is indicated by 
a filled circle. Dotted lines are loci of constant given 
by equation (||) ranging from lO^M© to lO^M©. The 
lines with their various core masses and maximum circular 
velocities are shown. They indicate that galaxies should 
rapidly loose all of their diffuse gas by ram pressure strip- 
ping. This figure shows that the gas in dwarf galaxies is 
completely stripped in rich clusters of galaxies within the 
short time-scale. 

The effect of radiative cooling has been neglected in 
our simulations, we estimate the role of radiative cool- 
ing in our models here. Since for low-mass dwarf galax- 
ies {Mq <^ lO^M©), the ram pressure already exceeds the 
gravitational force, radiative processes will not be able to 
aff'ect the gas stripping. However, it might play a role for 
massive dwarf galaxies {Mq ^ IO^Mq) which have deeper 
potential wells. The cooling time-scale is defined as (Efs- 
tathiou 1992) 



1 



fcRT 



"^cool — 



2^2(1-^)2 ngA' 



(35) 



where Y{~ 0.25) is the helium mass fraction, rig is the gas 
number density, and A is cooling rate. If we evaluate this 
timc-scalc for our standard model Mq — lO^M© using con- 
ditions immediately behind the leading shock, we find that 
Teooi = 3.18 X 10^ yr assuming T = 10^ K, A = 10-^3 erg 
cm~3^ and rig = 1.0 cm^^ from the result of our sim- 
ulations. This time-scale is larger than the instantaneous 
stripping time-scale (rjs = 2.03 x 10* yr) that is given by 
equation (p^). This implies that gas stripping from the 
massive dwarf galaxies is affected by the radiative cooling. 

However, we would then also need to consider the ef- 
fects of the star formation and subsequent feedback pro- 
cess such as stellar winds and supernovae heating from 



Mori & Burkert 



7 



massive stars. These feedback processes supply tlrernial 
energy to the gas and prevent efficient cooling (Dekel & 
Silk 1986; Mori at al. 1997; Mori et al. 1999). In ad- 
dition, the photoionization due to the ultraviolet back- 
ground radiation prevents gas cooling and keeps the gas 
hot (Efstathiou 1992; Babul & Rees 1992). Furthermore, 
though we studied only the interaction between the ICM 
and a hot interstellar medium in this paper, it is very in- 
teresting to examine also the fate of the cool and dense 
interstellar medium in a dwarf galaxy. In this case, the 
stripping history of the gas may be quite different even 
for the less-massive dwarf galaxies (M < IO^Mq). In a 
series of forthcoming studies, we will report the results of 
taking into account the multi-phase states of the gas with 
cool components in a galaxy, including the effect of the 
radiative cooling of the gas, energy input from stars, and 
the ultraviolet background radiation. 

We conclude that the ram pressure stripping of the dif- 
fuse gas from dwarf galaxies by the ICM is of prime impor- 
tance to understand the morphology-density relation for 
dwarf galaxies and the chemical evolution of dwarf galax- 
ies. Most of the heavy elements ejected by massive stars 
through Type II supernovae are found in the warm and hot 
diffuse gas. In order for these heavy elements to be incor- 
porated into stars, the hot gas has to cool and form cold 
clouds again. Our models indicate that this will be impos- 
sible for dwarf galaxies in clusters of galaxies, as the gas 
is stripped before condensing into clouds. Ram pressure 
stripping therefore plays a significant role for the chemical 
evolution of dwarf galaxies. Our results imply that most of 
the chemical evolution of dwarf galaxies must have occured 
at high redshift, before the virialized cluster had formed. 
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its hospitality throughout this research. Numerical cal- 
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ysis Center of the National Astronomical Observatory in 
Japan. 
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Fig. 1. — The density (upper) and the circular velocity (lower) of the Burkert's (1995) profile. Each curve is divided by the core mass of 
the dark matter halo which ranges from 10® Mq (left) to IO^^Mq (right). 

Fig. 2. — Snapshots for the model with core mass Mo = 10^ Mq, number density of the ICM ncc = 10~* cm~^, and relative velocity of 
the galaxy Ugai = 1000 km s"! as a function of elapsed time at 5.61 x 10®, 2.86 x 10^, 4.58 x 10^, 7.42 x 10^, and 9.63 x 10^ yr. The left and 
the right colour images show the logarithmic density distribution and the logarithmic pressure distribution, respectively. Arrows in the right 
panels indicate the velocity vector of the gas flow at each point. 

Fig. 3. — Snapshots for the model with core mass Mq = IQ-'^^Mq, number density of the ICM rtcG = 10^^ cm^^, and relative velocity of 
the galaxy v^^i = 1000 km s^^ as a function of elapsed time at 6.84 x 10^, 6.99 x 10*, 1.15 x 10^, 2.18 x 10^, and 3.26 X 10^ yr. The left and 
the right colour images show the logarithmic density distribution and the logarithmic pressure distribution, respectively. Arrows in the right 
panels indicate the velocity vector of the gas flow at each point. 

Fig. 4. — The evolution of the gas mass inside the core radius of the gravitational potential of the CDM halo as a function of time. 
Each column corresponds to a core mass of the CDM halo which varies from 10^ Mq to 1O^''M0. Each curve is categorized by the model 
parameters (dash-dot-dot-doted line: ncG = 10~^ cm~^ and v^^i = 1000 km s~^, doted line: ncc = 10~'^cm~^ and iig^i = 500 km s~^, 
dash-doted line: ncc = 10~^ cm~^ and Dgg,i = 1000 km s~^, dashed line: ncG = 10~^cm~^ and v^g^ = 500 km s~^, and solid line: 
ncG = 10~^ cm~* and Vg^i = 1000 km s~^). Each holizontal line indicates the time scale of lOrdyn. 

Fig. 5. — The mass loss time-scale for dwarf galaxies as a result of ram pressure stripping PCGi'gai due to the ICM as a sequence of the core 
mass Mo and the maximum circular velocity Uc.max of dwarf galaxies. Filled circles indicate the case that the gas is stripped within 1 Gyr. 
Open circles show that the gas is not stripped within 2 Gyr. Intermediate cases are marked by filed triangles. The shaded region corresponds to 
the inside of the median radius of the galaxy distribution for a typical galaxy cluster fitted by yS-model ncc = ncGo{l + (''Cg/^CGo)^}~^''^^i 
with /3 = 0.6, the central number density ncGO = 2 x 10~^ cm~^, the core radius rcGO = 0.25 Mpc, and the line-of-sight velocity dispersion 
(TQQ = 866 km s~^. tqq is the distance from the center of cluster. 

Fig. 6. — A scaled gas density is plotted against velocity dispersion for clusters of galaxies. Assuming a distribution of the gas that follows 
the /3-model, the gas number density is given at the median radius of the galaxy distribution. Data for the central gas density and /3 are from 
Jones & Forman (1999) and Brief, Henry &c Bohringer (1992), and the velocity dispersions are from these sources and from Hughes (1989). 
The filled circle corresponds to the Coma cluster. Solid lines are loci of constant Mct which is the critical core mass of the CDM halo, and the 
lines are labeled with their core mass and the maximum circular velocity. The lines indicate that dwarf galaxies should lose instantaneously 
all of their gas by ram pressure stripping. 
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